New interpretation for energy gap A of the cut-off approximation in the BCS theory 

of superconductivity 



O 

o 

(N 
C 



Zhidong Hao 

Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China 

(Dated: February 1, 2008) 

This paper concerns the solution of the self-consistency equation for energy gap parameter Ak in 
the BCS theory of superconductivity. We show that there exists a well-defined relation between the 
solution for energy gap parameter amplitude |Ak| for a general interaction Vk,k' and energy gap A 
obtained by using the cut-off approximation. The relation between |Ak| and A indicates that A is a 
weighted average over |Ak| of electronic states within cut-off energy £ c around the Fermi surface. In 
this interpretation for A, £ c is not a property of Vk,k', but a parameter specifying the energy range 
within which the weighted average over |Ak| is taken. We show that the proper choice for the value 
of £ c is only a few fesT c (i.e., £ c /kBT c is about 3 or 4). We also show that the cut-off approximation, 
even with £c/A;bT c = oo, is a good approximation when it is used to calculate quantities such as 
the condensation energy and the specific heat, but it leads to significant overestimation for the 
Josephson critical current density of a Josephson junction if £ c /kBT c S> 1 is assumed. 
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In the BCS theory of superconductivity,— i 2 - the super- 
conducting state is characterized by the existence of en- 
ergy gap par ameter Ak in quasi-particle excitation en- 
ergy i?k = \/£k + |Ak| 2 (where £k is the normal state 
electronic energy, measured relative to the Fermi level). 
Energy gap parameter Ak is determined self-consistently 
via the equation 



^ Vkk 2^ Ak ' 

k' K 



(1) 



where Vk,k' is the pairing interaction matrix element. 

In principle, once Ak is determined, thermodynamic 
quantities in the superconducting state can be quantita- 
tively calculated as functions of temperature T by start- 
ing with the diagonalized Hamiltoniar 1 ^ 1 ^ 



H = 



k 



7 kT 7kT + 7_ki7- 



kl 



(2) 



where U k = £k - £ k + |A k | 2 ta,nh(E k /2k B T)/2E k , and 
7ko- ancl 7ko- are the Fermi operators for quasi-particles 
in the superconducting state— i 2 - 

In general, Ak is a complex quantity, i.e., Ak = 
|Ak|e l£lk , and both amplitude Ak and phase #k can be 
wave vector k dependent. The cut-off approximation, 1 
in which Vk,k' is approximated by 



Vi 



k.k' 



-V < 







if |&|,|&' 
otherwise 



(3) 



and #k = is assumed to be a constant (which can be 
arbitrary), suppresses the k-dependence of Ak so that 



|A k | = 

and Eq. (Q]) becomes 



A > if |£ k | < £c 
otherwise 



1 



y £ 

|£k|<£c 



tanh ^k + A 2 /2fc B T 



2v^TA2- 



(4) 



(5) 



Cut-off energy £ c was thought to be of the same order 
as Debye energy Tllo d , i.e., £, c /k B T c ~ 7iw D /fc s r c > l, 1 
and ^c/ksTc — oo was often assumed in practical calcu- 
lation of various quantities—^ 

Major quantitative results of the BCS theory were first 
derived by using the cut-off approximation—^ Despite 
the fact that the approximation is oversimplified, the 
quantitative results have shown, in general, good agree- 
ment with experiments on a variety of (conventional) 
superconductors. 4,5 However, there are also noteworthy 
discrepancies. An example is that the predicted value 6 
for the magnitude of the Josephson critical current den- 
sity of a Josephson junction is much too large compared 
to what experimentally observed^u^ even though the 
prediction for the temperature dependence of the normal- 
ized Josephson critical current density has been found to 
be in excellent agreement with experiments. 

We have derived in Ref. [U a solution for energy gap 
parameter amplitude |Ak| for a general interaction Vk,k'- 
The solution for |Ak| shows that reduced energy gap pa- 
rameter amplitude |Ak|/fcsT c is a function only of re- 
duced variables ^kl/fcsTc and T/T c , which contains no 
explicit Vk,k'-dependence. The solution also shows that 
|Ak| is appreciable only for energies within a few ksT c 
around the Fermi level. This latter feature of |Ak| is 
very different from what one would expect from the cut- 
off approximation if ^c/fcsTe ^ 1 is assumed. Despite 
this difference, as we have shown in Ref. UH, the results 
for thermodynamic critical magnetic field H C (T), specific 
heat C (T) and normalized Josephson critical current den- 
sity I c (T)/I c (0), obtained by using the solution of Ref. 
[ill for |Ak|, are not much different from those obtained 
by using the cut-off approximation (with ^c/fceTc = oo). 
However, there is one significant difference: the value of 
2c(0), obtained by using the solution of Ref. [Ill for |Ak|, 
is only about a third of that obtained by using the cut-off 
approximation (with ^c/ksTc — oo). The reason behind 
these is further analyzed and made clear in this paper. 



2 



In the following, we show that there exists a well- 
defined relation between the solution for | | obtained in 
Ref. [Il| for a general interaction Vk,k' and energy gap A 
of the cut-off approximation. The relation between |Ak| 
and A indicates that A is a weighted average over |Ak| 
of electronic states within £ c around the Fermi surface. 
In this interpretation for A, cut-off energy £ c is not a 
property of the interaction, but a parameter specifying 
the energy range within which the weighted average over 
|Ak| is taken. We show that the proper choice for the 
value of £ c is only a few kgT c (i.e., £ c /£;_bT c is about 3 
or 4) . We also show that the cut-off approximation, even 
with ^c/ksTc = oo, is a good approximation when it is 
used to calculate quantities such as condensation energy 
H^(T)/8tt, specific heat C(T) and normalized Josephson 
critical current density I c (T)// c (0), but it leads to signif- 
icant overestimation for the magnitude of the Josephson 
critical current density if (.r./knTr 3> 1 is assumed. 

As we have shown in Refs. [U and [H, the following 
first order differential equation for |Ak(T)| holds: 



E^ 



dT 



tanh y/l- 



|Ak| 2 /2fc B T 



0. (6) 



This equation can be derived from Eq. ([TJ) by first op- 
erating d/dT on Eq. {!]), and then multiplying the re- 
sulting equation by (A£/2i?k) tanh(_Ek/2fcsT) and tak- 
ing summation over k. It can also be derived by calculat- 
ing entropy S from diagonalized Hamiltonian H of Eq. 
([2]), and letting the resulting expression for S to be the 
same as the standard expression expected for a system 
of Fermions i 11 ! 12 

Note that interaction Vk,k' and phase #k do not appear 
explicitly in Eq. Instead, critical temperature T c is 
involved through the condition that |Ak| = at T = T c . 
This indicates that |Ak| depends on Vk,k' only implicitly 
via T c . 

To see how critical temperature T c depends on inter- 
action Vk,k', we turn to Eq. (JTJ) . In the limit of T — > T c , 
we have |Ak| — * so that Eq. fT} can be linearized, and 
we have an eigenvalue problem: 



taah(|&|/2fc B r c 
2|&'| 



(7) 



In principle, critical temperature T c and phase #k are 
determined by solving this eigenvalue problem for given 
interaction Vk, k' and electronic energy spectrum £k- 

We turn back to Eq. ([6|) to consider how a solution for 
|Ak| can be obtained. Clearly, 



tanh (y^ + |Ak| 2 /2fc B T) tann( | a | /2 fc B T c ) 



16. 



(8) 



is a solution of Eq. ([5]). This equation, which was previ- 
ously obtained in Ref. |Tl|, is an implicit solution for |Ak| 



as a function of |£k| and T for given T c , and satisfies the 
condition that | A k | = at T = T c . 

However, Eq. JHJ is not the only possible solution of 
Eq. ©. Actually, as one can see, Eq. © can have in- 
finite number of solutions. For example, the solution of 
the form of Eq. ([4|) in the case of the cut-off approxima- 
tion is also a solution of Eq. ©. This can be seen by 
substituting Eq. into Eq. © to obtain 



E 



d 

dT 



tanh V£k + A 2 /2fc s T 



0, (9) 



and noticing that this equation can also be obtained from 
Eq. ([5]) by operating d/dT on it. Similarly, a solution 
of the form Ak = Awk, as in the case of a separable 
interaction Vk.k' = — VwkO'k'j is also a solution of Eq. 
([6]). We therefore need an additional constraint so that 
Ak(T)| can be uniquely determined. 

Note that diagonalized Hamiltonian H of Eq. @ is In- 
dependent because of its dependence on |Ak(T)|. This 
implies the existence of an additional self-consistency 
constraint, which, as we will see in the following, allows 
unique determination of |Ak(T)|. 

As we discussed in Ref. [TJ, since diagonalized Hamil- 
tonian H describes a set of independent quasi-particle ex- 
citations, there should be no coupling (except pair corre- 
lation) between the quasi-particle excitations. Therefore, 
we expect the thermal energy and the entropy associated 
with each pair of (kf, — k|) excitations to be 



£k = C/k + 2/ k £ k (10) 



and 



5 k = -2k B [/ k In / k + (1 - / k ) hi(l - / k )] , (11) 

respectively [where /k = (e Ek / kBT + l) -1 is the Fermi 
function]. However, as compared to these standard ex- 
pressions for £k and Sk, those derived from diagonalized 
Hamiltonian H of Eq. ^ contain additional terms in- 
volving dUk/dT, dEk/dT and df^/dT. Letting the sum 
of the additional terms in each expression to be zero, 
one gets a first order differential equation for |A k (T)|, of 
which the solution satisfying the condition |Ak(T c )| = 
is Eq. dHJiii Namely, the solution for |A k (T)| given by 
Eq. (|<3J) is the only solution that both satisfies the orig- 
inal self-consistency equation [Eq. (fT])] and ensures that 
Eqs. (fTt)|) and (fTTj) hold. [A complete solution for Ak is 
therefore a combination of the solutions of Eq. ([7]) for T c 
and 0k and the solution of Eq. (JHJ) for |Ak|.] 

There is a well defined relation between the quantity 
A of the cut-off approximation and |Ak| [hereafter, |Ak| 
means the solution for the energy gap parameter ampli- 
tude given by Eq. ©]. This can be seen as follows. 

Note that Eq. © is a first order differential equation 
for A(T), of which the solution satisfying the condition 
A(T C ) = is 



3 



E 

l?k|<«c 



tanh ( V^+A^/2fc B T) tanh (| & |/ 2 fc B T ) 



16. 



(12) 



By substituting Eq. |[5J into the above equation, we obtain 



E 

l?k|<Cc 



tanh y^k + A 2 /2fc B T tanh ^ + | A k | 2 /2fc B T 



= 0. 



(13) 



This equation defines A as a weighted average over |Ak| of electronic states with £k| < £ c ; we therefore can so 
interpret A. In this interpretation for A, £ c is not a property of the interaction, but a parameter specifying the energy 
range within which the weighted average over |A k | is taken. 

We next examine how A depends on £ c . From Eq. ([9]), by making the usual substitution $Z|£ k |<£ ~~ * N(0) ft dt; 
[where N(0) is the density of states at the Fermi level] and a rearrangement, we obtain the following expression: 



rf(A 2 
dT 



It « 


sech 2 


[y/e + A 2 /2T 






It 


tanh 


[ye + a 2 /2t) 


/ (e + A 2 ) 3/2 - sech 2 ( 


+ a 2 /2tJ 


/2T (£ 2 + A 2 ) 



(14) 



where energies are measured in units of k B T c . With this 
expression for d(A 2 )/dT and the initial value A 2 = 
at T = 1 (temperature T is measured in unit of T c ), 
we can numerically calculate A(T) for arbitrary £ c by 
using the Runge-Kutta method^ The integrals involved 
in the expression for d(A 2 )/dT are calculated by using 
the Simpson method.— 

We show in Fig. 1 the energy dependence of the gap 
parameter amplitude at T = 0. The results for A(0) 
obtained for £ c /k B T c = 1, 2, 3.37, 5, 10 and oo are shown 
as plots (a), (b), (c), (d), (e) and (f), respectively. (Here 
the number 3.37 is special, as we will see later.) Plot (g) 
shows |Ak(0)| versus |£k|- 

Also note that if a different cut-off range is assumed 
so that |Ak| = A for £ c i < |£k| < £c2 and zero otherwise, 
then A is simply a weighted average over |Ak| of elec- 
tronic states with £k in the range £ c i < |£ k | < £ C 2- As 
examples, the results for (£, c i/k B T c , = (6,7), 
(7, 8) and (8, 9) are shown as plots (h), (i) and (j), re- 
spectively, in Fig. 1. 

Comparing plots (a)-(f) and (h)-(j) with plot (g), it is 
evident that A is indeed a weighted average over |Ak| of 
energies within a specific range. This relation between A 
and |Ak| is mathematically expressed by Eq. (fT3"|) . We 
can also see from Fig. 1 that electronic states with lower 
|£ k | contribute with larger weights to the average. The 
fact that the difference between the values of A(0) for 
L/k B T c = 10 [plot (e)] and £ c /k B T c = 00 [plot (f)] is 
less than 1% indicates that contributions from electronic 
states with \£,\a\/k B T c 3> 1 are negligibly small. 

The minimum single quasi-particle excitation energy is 
-Ek,min = Ak F | (where kp is a Fermi wave vector). We 
have |Ak F (0)|/fc B T c = 2. The result A(0)/fc B T c = 1.764 
for £ c /k B T c = 00 [plot (f) in Fig. 1] is about 88% of 



-5c 




FIG. 1: Energy dependence of the energy gap parameter 
amplitude at T = 0. Plots (a)-(f): A(0)/fc s T c = 1.975, 
1.922, 1.859, 1.816, 1.778 and 1.764 for £ c /fc s T c = 1, 2, 
3.37, 5, 10 and 00, respectively. Plot (g): |Ak(0)| versus 
|f k |. Plots (h)-(j): A(0)/fc s T c = 0.517, 0.364 and 0.255 for 
{£,ci/k B T c ^ c2 /k B T c ) = (6,7), (7, 8) and (8, 9), respectively. 



\A kF (0)\/k B T c . 

In Fig. 2, we compare the T-dependence of nor- 
malized energy gap |Ak F (T)|/|Ak F (0)| with that of 
A(T) / A(0) . Four practically indistinguishable curves are 
shown in Fig. 2: |A k ,(T)|/|A k ,(0)|; and A(T)/A(0) for 
£,c/k B T c = 1, 100 and 00 of the cut-off approximation 
[A(T)/A(0) for t; c /k B T c = 00 was previously calculated 
by Miihlschlegel 3 ] . As shown in Fig. 2, A(T)/A(0) is 
practically ^-independent, and practically the same as 
|A kF (T)|/|A kF (0)|. 

We next examine how other quantities such as thermo- 
dynamic critical magnetic field H C (T), specific heat C(T) 
and Josephson critical current density I C (T) depend on 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

t = T/T c 



FIG. 2: Temperature dependence of the normalized energy 
gap. Four practically indistinguishable curves are plotted in 
this figure: |A kj? (t)|/|A kp (0)|; and A(*)/A(0) for £ c /k B T c = 
1, 100 and oo. 
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FIG. 3: Zero-temperature condensation energy _ff 2 (0)/87r 
versus cut-off energy £ c . The maximum _ff 2 (0)/87r = 1.613 
[in unites of N(0)(k B T c ) 2 ] is located at £ c /k B T c = 3.37; 
H 2 (0)/8iv — 1.556 for £ c /k B T c = oo. For comparison, the 
results calculated by using |Ak| is H 2 (Q)/8-k = 1.645. 

cut-off energy £ c when they are calculated by using the 
cut-off approximation. 

Once A(T) is obtained, quantities such as H C (T), C(T) 
and I C {T) can be calculated straightforwardly iia In Fig. 
3, we show the ^-dependence of zero-temperature con- 
densation energy H^(0)/8tt. The H 2 (0)/87r-versus-£ c 
curve shows a maximum iJ c 2 (0)/87r = 1.613 [in unites 
of N(0)(k B T c ) 2 } at ic/k B T c = 3.37. In the limit of 
£c/k B T c = oo, H 2 (0)/8n = 1.556. Note that these two 
values of H 2 (0)/8n are only about 2% and 5%, respec- 
tively, smaller than the result H 2 (0)/8ir = 1.645 calcu- 
lated by using Ak|. 

In Fig. 4, we show the T-dependence of thermody- 
namic critical magnetic field H c . The results are plotted 
as deviations from the 1 — (T/T c ) 2 law. For compari- 
son, the result calculated by using |Ak| and the results 
calculated by using the cut-off approximation for several 
different values of £ c are plotted in the figure. 

In Fig. 5, we show results for the T-dependence of 
electronic specific heat C. For comparison, the result cal- 
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FIG. 4: Deviation function D H = H C (T)/ H C {0)-[1-(T/T C ) 2 ] 
versus (T/T c ) 2 , calculated by using |Ak| and by using the 
cut-off approximation with ^different values of £ c /k B T c as in- 
dicated on the curves. 
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FIG. 5: Electronic specific heat C versus temperature T, cal- 
culated by using |Ak| and by using the cut-off approximation 
with different values of (, c /k B T c as indicated on the curves. 
The C(T c )/^/T c values for the four curves shown in the figure 
are 2.597, 2.525, 2.482 and 2.318, respectively. 



culated by using Ak and the results calculated by using 
the cut-off approximation for several different values of 
£ c are plotted in the figure as C/jT c versus T/T c [where 
7 = (27r 2 /3)fc|iV(0)]. At T = T c , we have C/jT c = 2.597 
by using |A k |, and C/-fT c = 2.525, 2.482 and 2.318 by 
using the cut-off approximation for ^ c jk B T c — 3.37, 100 
and 2, respectively. The result C(T c )/<yT c = 2.426 for 
£, c /k B T c = oo was previously obtained by Miihlschlegel,- 
From Figs. 3-5, we can see that, for calculating H c and 
C, the cut-off approximation gives results that are not 
much different from those obtained by using |Ak|. The 
approximation is optimized when £, c /k B T c = 3.37, sug- 
gesting the proper choice for £ c is about a few k B T c ^ 
We also note that the results for H c and C show only 
weak ^-dependence when £ c is a few k B T c or larger, so 
that the cut-off approximation remains a good approx- 
imation even with £, c /k B T c 3> 1. The reason for the 
weak ^-dependence of H c and C is that, as one can 
see from the expressions for H c and the relevant 

quantity for H c and C is quasi-particle excitation en- 
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FIG. 6: The ^-dependence of zero-temperature Josephson 
critical current density / c (0) of a symmetric Superconductor- 
Insulator-Superconductor junction. As indicated in the fig- 
ure, J c (0) = 0.834 and 2.771 (in units of k B T c /R n e) for 
S^c/ksTc — 3.37 and oo, respectively. The result obtained 
by using |A k | is J c (0) = 1.103. 
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FIG. 7: Comparison between the T-dependence of nor- 
malized Josephson critical current density 7 C (T)/7 C (0) of 
a symmetric Superconductor-Insulator-Superconductor junc- 
tion calculated by using |Ak| and those by using the cut-off 
approximation with different values of ^ c /fcsTc as indicated 
in the figure. 

ergy T k = ^ + A 2 (or E k = ^ + |A k | 2 ), which 
becomes Ek — |£k| = independent of A (or |Ak|) for 
|£k| 3> ksT c , so that only electronic states within a few 
UbT c around the Fermi surface contribute significantly to 
the difference between the superconducting and normal 
states. The situation is different in the case of Joseph- 
son critical current density I c of a Josephson junction, 
because, as one can see from the expression for I c r^ the 
relevant quantity for I c is energy gap parameter Ak it- 
self. This leads to a strong ^-dependence for I c when it 
is calculated by using the cut-off approximation, as we 
will see next. 

In Fig. 6, we show the ^-dependence of zero- 



temperature Josephson critical current density T(0) of 
a symmetric Superconductor-Insulator-Superconductor 
junction. As shown in the figure, T(0) is a monoton- 
ically increasing function of £ c that does not saturate 
until £ c //cbT c ~ 10 3 . As indicated in the figure, we 
have 4(0) = 0.834 and 2.771 (in units of k B T c /R n e) for 
^c/ksTc = 3.37 and oo, respectively (the result of I c for 
^c/ksTc = oo was previously obtained by Ambegaokar 
and Baratoff 6 ). For comparison, the result T(0) = 1.103 
calculated by using |Ak| is also indicated in the figure, 
which is only about 40% of the result calculated by using 
the cut-off approximation with ^ c //csT c = oo. It is evi- 
dent from the results shown in Fig. 6 that the cut-off ap- 
proximation significantly overestimates I c if £ c /fc_BT c 1 
is assumed. 

In Fig. 7, we show the T-dependence of normalized 
Josephson critical current density I c / I C (Q) of a symmetric 
Superconductor-Insulator-Superconductor junction. For 
comparison, results calculated by using the cut-off ap- 
proximation for several different values of £ c and that by 
using |Ak| are plotted in the figure. As we can see from 
the figure, for lower values of £ c /&bT c , the difference be- 
tween the result obtained by using the cut-off approxima- 
tion and that by using |Ak| is appreciable. However, the 
difference becomes much less significant for £ c /A;bT c 1. 

Experimentally,— ^ i 10 ' 16 i 17 the T-dependence of 
J c /J c (0) obtained by using the cut-off approximation 
with S^c/ksTc = oo, which is practically indistinguishable 
from or very close to that obtained by using |Ak|, has 
been well confirmed, but the magnitude of I c has 
always been found to be much too small compared to 
that predicted by using the cut-off approximation with 
^c/ZcbT; = oo. Clearly, as we can see from Figs. 6 and 
7, the results obtained by using |Ak| can help improve 
agreement between theory and experiments. 

In summary, we have shown in this paper that energy 
gap A that was originally obtained by using the cut-off 
approximation can be interpreted as a weighted average 
over |Ak| of electronic states within £ c around the Fermi 
surface [here |Ak| is the energy gap parameter amplitude 
for a general interaction Vk.k', and satisfies Eq. ©]. In 
this interpretation for A, cut-off energy £ c is not a prop- 
erty of the interaction, but a parameter specifying the en- 
ergy range within which the weighted average over |Ak| 
is taken. We have shown that the proper choice for the 
value of £ c is only a few ksT c (i.e., ^ c /fcsT c is about 3 
or 4). We have also shown that the cut-off approxima- 
tion, even with £ c /fc_BT c = oo, is a good approximation 
when it is used to calculate quantities such as condensa- 
tion energy H^(T)/8tt, specific heat C(T) and normal- 
ized Josephson critical current density 7 C (T)// C (0), but 
it leads to significant overestimation for the magnitude 
of the Josephson critical current density if ^ c /^sT c ^> 1 
is assumed. 



1 J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 2 N. N. Bogoliubov, Nuovo Cimento 7, 794 (1958); Zh. Eksp. 
108, 1175 (1957). 



6 



Teor. Fiz. 34, 58 (1958) [Sov. Phys. JETP 7, 41 (1958)]; 
J. G. Valatin, Nuovo Cimento 7, 843 (1958). 

3 B. Miihlschlegel, Z. Phys. 155, 313 (1959). [English trans- 
lation in The Theory of Superconductivity, edited by N. N. 
Bogoliubov. Gordon and Breach, Science Publishers, New 
York (1968).] 

4 J. Bardeen and J. R. Schrieffer, in Progress in Low Tem- 
perature Physics (Edited by C. J. Gorter) Vol. 3, p. 170. 
North-Holland, Amsterdam (1961). 

5 R. Meservey and B. B. Schwartz, in Superconductivity 
(Edited by R. D. Parks). Marcel Dekker, New York (1969). 

6 V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 
(1963); erratum 11, 104 (1963). 

7 P. W. Anderson and J. M. Rowell, Phys. Rev. Lett. 10, 
230 (1963). 

8 M. D. Fiske, Rev. Mod. Phys. 36, 211 (1964). 

9 I. K. Yanson, V. M. Svistunov, and I. M. Dmitrenko, Z. 
Eksp. Teor. Fiz. 47, 2091 (1964) [Sov. Phys. JETP 20, 
1404 (1965)]. 

B. D. Josephson, in Superconductivity (Edited by R. D. 
Parks). Marcel Dekker, New York (1969). 



11 Z. Hao, Mod. Phys. Lett. B7, 1439 (1993). 

12 Z. Hao, J. Phys. Chem. Solids 57, 1215 (1996). 

13 S. D. Conte and C. de Boor, Elementary Numerical Anal- 
ysis: An Algorithmic Approach, 3rd edition, McGraw-Hill, 
New York (1980). 

14 See Ref. or references cited therein for details about ex- 
pressions for H c , C and I c . In calculating these quantities, 
as usual, the substitution ^ k — > ^(0) J d£ is made, and 
resulting integrals are numerically calculated by using the 
Simpson method.— 

The value ( c /feT c = 3.37 is obtained by maximizing con- 
densation energy H^/8tt at T = 0. The optimize value for 
£ c may show some degree of T-dependence, but it is not 
important, as our point here is to show that the proper 
choice for £ c is only a few fcsT c . 

16 E. P. Balsamo, G. Paterno, A. Barone, P. Rissman, and 
M. Russo, Phys. Rev. B 10, 1881 (1974). 

17 A. G. Sun, D. A. Gajewski, M. B. Maple, and R. C. Dynes, 
Phys. Rev. Lett. 72, 2267 (1994). 



